Feedback dynamic control for exiting a debt-induced spiral in a deterministic Keen model

The Keen model is designed to represent an economy as a dynamic system governed by the interactions between private debt, wage share, and employment rate. When certain conditions are met, the model can lead to a debt spiral, which accurately mimics the impact of a financial crisis on an economy. This manuscript presents a recipe for breaking this spiral by expressing Keen’s model as an affine nonlinear system that can be modified through policy interventions. We begin by considering critical initial conditions that resemble a financial crisis to achieve this goal. We then locate a desired point within the system’s vector field that leads to a desirable equilibrium and design a path towards it. This path is later followed using one-step-ahead optimal control. We illustrate our approach by presenting simulated control scenarios.


Introduction
The global financial crisis of 2008 stands as a significant event in economic history.It was sparked by the collapse of the US subprime mortgage market, leading to severe repercussions across the world economy.Researchers have extensively analyzed this and other crises in financial markets [1][2][3], providing valuable insights into their causes and consequences.According to studies by Reinhart and Rogoff [4], one of the underlying reasons of the crisis was the rapid accumulation of excessive household debt, fueled by lax lending standards and financial deregulation.Minsky [5,6], had previously warned about the inherent instability of financial systems and the tendency for economic cycles to swing from periods of exuberance to crises [7][8][9].This concept, known as the "financial instability hypothesis" (FIH), resonates with the dynamics observed during the crisis, which ultimately led the economy to a debt spiral with high unemployment and low wages.
Re-directing the economy to a stable equilibrium can be seen as an engineering control problem in a dynamical system [10][11][12].The specialized literature [13][14][15][16][17] has shown the importance of incorporating control theory in economic and financial analysis, highlighting its potential to enhance policy effectiveness and mitigate economic fluctuations [18,19].Furthermore, the application of control techniques from different fields [20][21][22][23][24] has a long tradition in macroeconomics [25,26], where models incorporating feedback control [27] have been employed to address issues such as inflation control [28], and monetary policy implementation [29].Moreover, control techniques have found utility in the study of chaotic financial systems [30,31], and in applications for risk assessment, portfolio optimization, and asset pricing, see e.g.[32].For example, research by Zheng et al. [33] has demonstrated the benefits of using control theory to develop robust investment strategies that adapt to changing market conditions and enhance financial performance; while the works of [34,35] have shown the use of control engineering to manage private and public debt respectively.
In this document, we follow Minsky's FIH and use Keen's model [36,37] to characterize the economic system.This model emphasizes the interplay between debt, asset prices, and firm behavior that occurs during a financial crisis, mimicking the non-linear cyclical behavior of the economy.Our main objective is to apply engineering control techniques to the Keen model and present a recipe to break the debt spiral by expressing it as an affine nonlinear system that can be modified through policy interventions.We begin by considering critical initial conditions that resemble a financial crisis to achieve this goal.We then locate a desired point within the system's vector field that leads to a desirable equilibrium and design a path towards it.This path is later followed using one-step-ahead optimal control, which situates this study in the realm of dynamically controlled systems [38,39].
We have formalized our findings in the form of propositions and theorems to summarize our ideas and justify our proposals.Our main result (see theorem 3.1 in section 3) shows that there exists a unique control law for the Keen model.This result allows us to present a roadmap to exit the debt spiral characterizing the financial crisis of 2008.We stress that the utilized control policies require the crucial assumption that certain parameters in the Keen model can be altered or "controlled" under a certain level of latency.Namely, we consider that it is possible to control real interest rates with almost no latency, as its level is a decision made by the Central Bank; but adding extra dynamics to allow for a slow response when controlling the growth rate of labor productivity and the growth rate of the labor force, which, if controllable, can only be assumed to change at a much slower rate.To the best of our knowledge, this represents the pioneering integration of control engineering techniques into the Keen model.Consequently, this research augments the theoretical and empirical discourse surrounding this specific financial system in academic literature.
The rest of the manuscript is organized as follows.In section 2, the Keen model is presented.In section 3 the Keen model is revisited and is presented in the form of a non-linear affine system and the conditions of its controllability are discussed.In section 4, a simulation study is performed to explore the behavior of the system under different stress scenarios.Finally, Section 5 provides a discussion of the results and the conclusions resulting from this research.To make this a self-contained document, we added a technical appendix summarizing the basic ideas of linear control systems and feedback linearization that utilized in this manuscript.

Modeling financial instability
A state-space control model represents a dynamical system in terms of its states, inputs, and outputs.In economic terminology, the states are known as endogenous variables, while the inputs are referred to as exogenous variables.The outputs, on the other hand, do not have a specific name within economics, but they can be considered as performance variables that relate to a predetermined target for the system.For example, in the context of the money market, the Central Bank aims to decrease inflation (output) by increasing endogenous interest rates (state) via the reduction of exogenous money supply (input).
In this section, we present the Keen model [36] as an appropriate representation of the economy that enables us to demonstrate a debt share control mechanism in the system's state space.We use this terminology to describe a procedure for guiding the current state of the economy from an initially unstable region to a desired location in a stable region.While exact tracking of the reference trajectory is not our primary concern, it is essential that the system states are not significantly distant from such a reference and do not violate the positivity of the interest rate.

Goodwin Model
Consider the Goodwin model [40] as a representation of an economy, resulting from the interaction between the wage share (x 1 ) and the employment rate (x 2 ), governed by the following system: where α is the rate of growth of labor productivity, β is the rate of growth of labor force, λ is the capital depreciation constant, ν represents the capital-to-output ratio, and FðxÞ ≔ � 1 ð1À xÞ 2 À � 0 , denotes the Phillips curve for constants ϕ 0 , ϕ 1 > 0. In addition, note that the state variables x = (x 1 , x 2 ) > take values in the set (0, 1) × (0, 1), where (0, 1) denotes the open interval between 0 and 1. Enforcing the condition it is clear that (1) has a unique non-hyperbolic stable equilibrium and its solutions are given by periodic curves centered at (3).Fig 1 shows the phase portrait of the system around (3) and illustrates the endogenous cycle of the system.We note that the Goodwin model also allows for an unstable equilibrium at (0, 0) and that this point is a saddle point if: i) F(0) − α > 0 and ii) 1 n À a À b À l < 0. However, we do not study this equilibrium further as it is of no economic interest.
Model (1) can also be studied from a control perspective.We call a system controllable if there exists a policy, in the form of inputs to the dynamical system, that can move it from an initial state to a final one, see e.g.[38,39].Here, we will briefly assess the controllability of the Goodwin model for inputs u 1 ≔ α and u 2 ≔ β, in terms of the system's Taylor linearization and a weaker form of controllability known as accessibility (see Appendix A).We stress that our intention is to provide the necessary arguments to illustrate that (1) can indeed be manipulated, see e.g., [41,42], and not to present a complete treatment of the theory of controllability of dynamical systems.
For a Taylor linearization procedure, let x 0 = (x 01 , x 02 ) > and u 0 = (u 01 , u 02 ) > give nominal points around which a Taylor expansion of (1) can be performed, such that x � x 0 + δx and u � u 0 + δu.One thus has the linearized dynamics where Therefore, from (27), one has that where it is clear that RankðC A;B Þ ¼ 2 for any x 0 2 (0, 1) × (0, 1), which ensures the existence of  1.The first quadrant is divided into four regions according to the growth and acceleration of wage share (x 1 ) and employment rate (x 2 ), around the stable non-trivial equilibrium given at (0.835, 0.969).(b) Simulated trajectories of wage share (x 1 ) and employment rate (x 2 ) for the system in (a). https://doi.org/10.1371/journal.pone.0295859.g001 a local input policy u = (u 1 , u 2 ) > that can manipulate (1) towards a desired state.More generally, one can check the accessibility of (1) (see Appendix A), which we present in the form of the next proposition.Proposition 2.1.The Goodwin model described in (1) is accessible for all states x in the set (0, 1) × (0, 1).
Proof.Showing that (3) is accessible follows from computing (32) and checking that its rank is 2. That is, from (1), one can see that 5 where it is clear that Rank(Q(f, g)) = 2 for all x 2 (0, 1) × (0, 1).Hence the system is locally accessible in (0, 1) × (0, 1), which is the region of interest for the control of the Goodwin model, and thus it is accessible.

Keen model
The Keen model can be seen as a direct extension of the Goodwin model in the following differential system.
where x 3 represents the economy's debt share, r is the debt's real interest rate, and κ is a private investment function.Note that the first two equations in (4) are identical to those in (1) for the appropriate choice of κ(�).In (4), we define this function as κ(x 1 , x 3 ) ≔ κ 0 + κ 1 tan −1 (κ 2 (1 − x 1 − rx 3 ) + κ 3 ), for constants κ 0 , κ 1 , κ 2 , κ 3 > 0. Note that the investment function increases with respect to the net profit, i.e. profit after wages and debt are paid It is easy to see that model (4) has a unique equilibrium given by where � z ¼ 1 À � x 1 À r� x 3 .Moreover, one can show that the equilibrium in ( 5) is stable as long as � x 1 > nl, but might not be finite.This is the main point of the Keen model, which given certain initial conditions, allows for a stable equilibrium with infinite debt, zero wage share, and zero employment rate.Hereafter we denote R u the set of initial conditions in R 3 þ that lead to this undesired equilibrium, and R d to its complement, thus related as Similar to the procedure followed for the Goodwin model, one can select α, β, and r as the inputs to which policymakers can assign values for altering the dynamics of (4).However, the system will not have affine form as in ( 29), which will be necessary for constructing an appropriate control policy.Therefore, we create a new state x 4 ≔ r, and select u 1 ≔ α, u 2 ≔ β and u 3 ≔_ r as the inputs that policymakers can modify at will to alter the dynamics of the system.The new state satisfies which implies that one changes the variation of the interest rate rather r itself.This augmentation of (4) has the form of a fourth-dimensional control affine system as in (29).The next proposition shows that the Keen model is accessible.
Proof.Showing accessibility is accomplished by computing (32) and checking that its rank is 4.That is, one can see that ; where þ , which makes the augmented Keen model accessible in the domain of interest.
For clarity of the exposition, all parameter values used in this manuscript for the Goodwin and the Keen models are listed in Table 1.: We close this section stressing that the adequacy of both the Goodwin and Keen systems to model economic phenomena have been extensibly explored.We refer the interested reader to [43][44][45][46][47] for the Goodwin model, and to [48,49] for the Keen model.

Dynamic control to scape a debt spiral
Our objective in applying control theory to the Keen model is to steer the system from x 0 2 R u to any point in the stable region x τ 2 R d , as any point in this region will take the system to the desired finite equilibrium.See Table 2 for the exact values of these points.
This goal can be achieved by tracking some oracle reference path ℓ with initial point ℓ(0) = x 0 and final point ℓ(τ) = x τ with the aid of a given control mechanism, called control policy.This controlled path can be seen as the solution to a problem that minimizes the gap between the current state of the system and the oracle path at each step in time.

Input variables
Input variables or control policies, denoted as vector u, serve as the means of interaction between a policy maker and the economic system described in (4).For the Keen model, we have selected inputs indirectly related to parameters α, β, and r, in the following manner: where r 0 = 0 and α 0 and β 0 take the values used for the illustration of model ( 4) in Table 1.In addition, the values τ 1 , τ 2 > 0 are set to calibrate their variation latency.Note that input u = (u 1 , u 2 , u 3 ) > interacts with the system trough new states x 4 , x 5 and x 6 .In the case of x 4 we consider an immediate control over the change in the real interest rate, but for the cases of α and β we consider a delayed effect as described by ( 7) and ( 8) respectively.Hence, system (4) can be re-written isolating the effect of u over the system as an affine nonlinear system given by: which can be explicitly written as   where κ(x) has the same meaning as before, but given that x 4 = r in ( 6), we can re-write it as κ . System (9) makes a clear distinction between input variables u, state variables x, and output variables y.For the Keen model, this means that the system can be controlled by modifying inputs u given by the growth rates of labor productivity (x 5 ), labor force (x 6 ), and interest rates (x 4 ).The impact on these quantities is then carried into the system's full state x, out of which only certain output quantities y, given by the state variables in

Exact linearization
To control the inputs in model (9) it is convenient to first linearize it.Here we employ a technique called feedback linearization, which allows for an exact linearization of the system, in contrast to e.g.Taylor linearization, which is only an approximation.The main idea of the method is to find a feedback linearizing control law (see Eq (42) in subsection B.2 of the appendix) that linearizes an affine nonlinear system exactly.The following theorem shows that for the particular case of system (9) this control law exists.
Theorem 3.1.There exists a unique feedback linearizing control law u K (x) for system (9) given by where In addition, v and Θ K (x) are three dimensional vectors, where the first is a free vector and the second has elements given by where we have omitted subscript K for each component to ease the notation.More specifically À 0:16Þ À 0:5Þ; Proof.The proof follows the construction of (42) as provided in subsection B.2 of the Appendix.The computations of D K (x) and Θ K (x) follow from ( 39) and ( 41) respectively.To make (9) exactly linearizable one must show that D K (x) is invertible in the domain x 1 2 (0, 1), x 2 2 (0, 1), x 3 2 R and x 4 2 R such that (12) holds, which can be trivially achieved given that x 4 is controlled by the modeler via u 3 .

Controlling the Keen model
Consider an Euler discretization of (14) for some fixed sampling time δ as where M, N, and W where defined in (15), and d 2 R represents some sampling time.

Given certain reference signal ℓ[k], one can define a quadratic loss function for error e[k] = y[k] − ℓ[k]
as where Q is a weight matrix affecting the coordinates of the error.Note that function ( 17) can be expressed explicitly in terms of the control policy u[k] by replacing vector v[k] from (10) into the definition of y[k] in (16).Thus, we can re-write it as where the loss function is conveniently expressed in terms of the control policy vector.In addition, to bound the search space of the policy vector, we consider the following constraints: i) |u i | � B u , i = 1, 2, 3 where B u is some real scalar; and ii) r min < r < r max , which can be written in terms of u i as (r min − x 4 )/δ � u 3 � (r max − x 4 )/δ, where δ.With the latter considerations, we can express our final optimization problem as which is convex in u[k], and thus existence and uniqueness of the solution are guaranteed.

Simulation results
This section presents a simulation study to showcase the methodology presented in the previous section.This synthetic experiment utilizes the parameters specified in Table 2, following [36], thus the use of real data remains open for future research.We begin by providing an algorithm to generate a reasonable tracking trajectory for the Keen model ( 4); continue presenting two simulation scenarios, and close by inspecting some numerical results.

Reference path
To obtain an appropriate reference path, a gradient-based planner (GBP) algorithm is introduced, which explicitly embeds Keen's model dynamical equations in the path generation between two pre-established locations in the state space.Let x 0 2 R u and x τ 2 R d , where R d and R u denote desired and undesired regions in the state-space respectively.Let F x denote the vector field of (4) and F K denote an attractive potential field with attractive constant K centered at x 0 , that strengthens with the proximity to x τ .Define the reference path as the minimum resistance path from x 0 to x τ under the force field generated by The construction of such reference path can be recursively performed as presented in Algorithm (1).
Algorithm 1 Gradient Based Planner It is worth noting that there is a trade-off in the construction process of the reference path.Namely, while a shorter path might be desirable, it does require more control effort.This trade-off should can be managed by the selection of K, in a way that the path is short but the effort is feasible.Another important factor is the singularity of the system at x 2 = 1 due to the Phillips curve.If the attractive force is not strong enough, the path generated by F s could enter this singularity and never leave.In this study, the value of K was selected balancing control effort and length of the reference path.

Simulation scenarios
The control policies enter the system in the form of u 1 , u 2 , and u 3 , representing changes in the rates of growth of labor productivity (α), labor force (β) and interest rates (r) respectively, starting from a certain initial state of the economy x 0 2 R u .From this point onwards the control policies presented in section 3 are utilized to transfer the economy from x 0 towards R d .We selected the specific point x 0 = (0.5, 0.5, 0.8) > considering that, while in the undesired region, it remains close to the frontier with the desired region, which is easy to note by inspecting Fig 3, where the cases for x 3 = 0.1 and x 3 = 5 are represented in panels (b) and (c) respectively.
The reference path in Fig 4 is generated utilizing Algorithm 1 in Subsection 4.1, so that control can be applied to track it.The minimization procedure described in Subsection 3.3 minimizes the error between output and a reference path using the following parameter values: B u = 1, r min = 0, r max = 0.5.Furthermore, we considered two main illustrative scenarios.In the first case, we presumed low latency (τ 1 = τ 2 = 1000) in any of the controlled variables, and in the second case, we assume high latency (τ 1 = τ 2 = 1), with the exception of interest rates.We stress that the second case is more realistic as control policies u 1 and u 2 are not expected to be executed with immediacy, and include the first scenario for comparison purposes exclusively.

Results
We illustrate the control process in  In fact, the reference path (continuous black line) and the controlled path (gray dashed line) are very close to one another, indicating that the system can be controlled with ease.In scenario (b) we observe that the final (black) and target (gray) points are further from each other and that the controlled path does not follow the reference path with the same proximity as in the previous scenario, which resulted in a slower response to the applied control policies due to the added latency.Nonetheless, in both scenarios, we see that the Keen model is effectively controlled from an initial point (x 0 ) towards a small enough neighborhood of final destination (x τ 0 ) within the stable region of the Keen system.Once the control policy is released in such neighborhood the uncontrolled path follows a stable downward (continuous gray line). https://doi.org/10.1371/journal.pone.0295859.g007 account for the fact that the target might only be reached within a certain neighborhood, which is good enough for our purposes as long as x 0 t 2 R d , which is indeed the case.Fig 8 shows the trajectories of the resulting levels of α, β, and r for the two latency cases in panels (a) and (b).The main difference between the two cases can be seen in the behavior of interest rates.While in the heavy latency case, the goal can only be achieved by increasing interest rates to levels surpassing 0.4, in the low latency case, interest rates are modified within a narrower radius.Besides the increment of interest rates, the recipe for leaving a debt spiral seems to pass through decreasing the rate of growth of labor productivity and increasing the rate of growth of the labor force, which might seem counter-intuitive at first sight.However, since Keen's model inherits Goodwin's trade-off behavior between firms and households, the real meaning of this behavior is that either production must be reduced, and more of the utility should go to households or more labor should be hired.
We stress that the success of the former control strategy assumes that the Keen model provides a complete characterization of the economy, and it is conditioned to the availability of control policies in interest rates (u 3 ) and the growth rate of labor productivity (via u 1 ), and the growth rate of the labor force (via u 2 ).Lastly, for the former recipes to be successful, they must be deployed exactly as described, starting from point x 0 = (0.5, 0.5, 0.8) > , under the aforementioned latency conditions.

Conclusions
This manuscript shows that the Keen model can be controlled to steer the economy from an initial undesirable location, where debt grows very fast, to a desired location where all states converge over time to a finite, and well-balanced equilibrium.To accomplish this goal, the model was dynamically extended and feedback linearized so that one-step-ahead optimal control can be applied.It is thus argued that the linearized version of the Keen model can be used as a baseline reference to obtain policies that control debt share within an economy in critical  and (b) shows the same for the case of high latency (τ 1 = τ 2 = 1).As can be seen, in the low latency case, the control policies are able to fluctuate more freely than in the high latency case, as the effect of the latter is the reduction in the velocity of the policies.In both cases, however, is noticeable the reduction in the rate of growth in labor productivity, as well as the increment in the rate of growth of the labor force, while maintaining positive values for the interest rates, as enforced by the conditions in equation (20).https://doi.org/10.1371/journal.pone.0295859.g008situations.We refer to our findings as a baseline reference as the success of the suggested policies depends crucially on the mild controllability of the productivity rate of growth, and the labor force rate of growth, which are assumed.
It is worthwhile noting that under the Keen model, a financial crisis is not considered an anomaly, or a black swan, but a cyclical non-periodic phenomenon, that will, unequivocally, repeat itself.In this sense we consider our results to be of great significance for high-level policy-making since they allow for a deeper understanding of the financial system and the control policies that can be utilized to navigate the next financial crisis.In particular, theorem 3.1 shows that there exists a unique control law for the Keen model, which can be used to exit a crisis under certain specific conditions.This work can be extended in various directions.From a statistics perspective it is worth noting that the empirical testing of the Keen model has not been as vastly evaluated as in the case of the Goodwin model, and systematic validation of the Keen model is foundational for any future developments.From a more engineering perspective, the reference generation algorithm can be improved by considering control action and effort as it is done in path planning.Furthermore, control effort can also be included in the cost function of the minimization problem for the policy generation, so that input variability can also be controlled.Lastly, the stochastic and data-driven extensions are currently under investigation and are deferred to future publications.

Appendix
This appendix presents the basic ideas of linear control systems and feedback linearization that are required to follow the steps in section 3. The contents are based on [21,38,39], and include linear and geometric-based techniques.
In general, a control system with states x 2 R n and inputs u 2 R m is any system that has the form In addition, the system is called controllable if for any given x 0 ; x t 2 R n and t 2 R, there exists an input u : R !R m such that the problem has a solution.

A Control systems and their controllability
A control system is called linear control system if it can be written as where x 2 R n , u 2 R m , y 2 R q and A, B and C are matrices of appropriate dimensions.
Additionally, we say that the system is controllable if the rank of its controllability matrix equals n.This implies that n directional vectors can be employed to steer x to any desired place in R n .A control policy u = Kx, K 2 R m�n is known as state feedback, and it provides the desired policy by selecting the components of K appropriately.If such K exists, the system (26) turns into and its trivial solution is x = C exp{(A + BK)t}x 0 , which can be made zero by choosing K in a way that the eigenvalues of A + BK are negative and large enough.Hence, the system is guaranteed to reach the origin after some time commanded by how fast the exponential decays, which is a tuning parameter for the control policy.When a control system is not linear, it is still possible to Taylor linearize it so that one can employ techniques developed for linear systems, however, such Taylor linearization is only valid over a small neighborhood of an operating point.This limits the applicability of the control policy, the quantification of uncertainty is not accurate, and no guarantees of performance can be claimed.
In a more general setting, hereafter we consider the class of affine nonlinear systems.That is, systems of the form with x 2 R n and f, g and h sufficiently smooth in a domain of R � R n .

A.1 Controllability of affine nonlinear systems
Bringing in concepts from differential geometry [50], the set of all smooth vector fields on a state space is denoted by XðRÞ.This forms a real vector field and a Lie algebra under the multiplication defined by the bracket operation ½�; �� : XðRÞ � XðRÞ !XðRÞ For each constant input, a vector field is defined by (21), and the collection of such vector fields is denoted by V 0 ðRÞ.The smallest subalgebra of XðRÞ containing V 0 ðRÞ is denoted by VðRÞ.Furthermore, denote by TVðRÞ x the space of tangent vectors spanned by the vector fields in VðRÞ at x.A system is said to satisfy the controllability rank condition at x 0 if the dimension of TVðRÞ x 0 is n.
However, not all vector fields are relevant to construct TVðRÞ x .Here it will be enough to focus on vector fields built iteratively by brackets of the form where ad 0 f ¼ f .Thus, for an affine system of the form _ x ¼ f ðxÞ þ P m i¼1 g i ðxÞu i the controllability rank condition at x reduces to checking the rank of the matrix equals n for some m > 0. When this condition is satisfied, then such a system is called accessible.The property of accessibility simply means that the set of states that can be reached from a state x through some admissible input is not empty.Therefore, stronger notions of controllability exist for nonlinear systems, but these are outside the scope of this manuscript.More details and stronger notions of controllability for nonlinear systems can be found in [51][52][53][54].Nevertheless, one can see that (32) reduces to (27) for linear systems as in (26) due to the Cayley-Hamilton Theorem [39].
In the following section, we present a technique that guarantees the manipulation of systems of the form (29) avoiding Taylor linearization, called feedback linearization.This method allows for the design of a matrix K, as in (28), for tracking a desired reference path without resorting to approximations.

B Feedback linearization
Feedback linearization is a technique used to obtain a linear system that exactly represents the dynamics of an affine system as presented in (29).To accomplish this goal, a nonlinear change of variables is needed.The procedure to exactly linearize an affine system is described next.Consider first the definition of a Lie derivative.
Definition 7.1.Let f and h be smooth functions.Denote by L f hðxÞ ¼ @h @x f ðxÞ the vector field Lie derivative of h with respect to f, and is the k-th iterated Lie derivative of h with respect to f.The relative degree of an affine nonlinear system is then provided by the next definition.Definition B.2. System (29) is said to have relative degree 0 � p � n in the region R x 0 2 R n if for all x 2 R x 0 L g L iÀ 1 f hðxÞ ¼ 0 for i ¼ 1; 2; . . .; p À 1;

B.1 SISO systems
For a system with a single input and a single output (SISO), the k-th derivative of y 2 R in (29) for a neighborhood of t 0 = x(0) can be written as and for relative degree p > k it holds that y ðkÞ ¼ L k f hðxÞ.Therefore, the relative degree p of a system is equal to the number of times the output must be differentiated to have the input After this change of variables, the resulting system is: We say that an affine system is exactly linearizable if there exists an invertible

u: ð37aÞ
In this case, the system has the form (26), and its controllability matrix is given by which has full rank and thus the system is controllable.We highlight the generality of this result, which means that if the original system is exactly linearizable (D is invertible), then the resulting linear system is always controllable.

B.2 MIMO systems
For a system with ℓ outputs and m inputs, one can compute the relative degree for each output, say p i for i = 1, . .., ℓ, by simply computing derivatives until any of the inputs appear.This defines the vector relative degree (p 1 , . .., p ℓ ), which, once again, is well-defined if is invertible.If such is the case, then D is known as the system's decoupling matrix.Even in the case of a non-square matrix which means that the number of outputs and inputs are not equal (ℓ 6 ¼ m), the pseudo inverse can be taken [55].Analogous to the case of only one output, this creates ℓ decoupled linear systems in normal form.Let z i,j be the j-th component of the i-th system, the linearized system in Byrnes-Isidori normal form is given by for i = 1, . .., ℓ.As before, this system was obtained from (33), where one can define a chain of integrators as above.Defining and putting together (33) for all outputs gives the control law Note that ( 40) is also a linear system that is controllable.A standard state feedback law can now be applied as in (28) so that the overall system is stable and can be made to follow a desired path.The state feedback law is applied on v i = K i z i , where z i ¼ z 1;i � � � z 1;p i Þ > À can be obtained from the derivatives of each output and K i is a user-defined vector of gains.

Fig 1 .
Fig 1. Goodwin model.(a) Phase portrait trajectory of the Goodwin model for the parameters in Table1.The first quadrant is divided into four regions according to the growth and acceleration of wage share (x 1 ) and employment rate (x 2 ), around the stable non-trivial equilibrium given at (0.835, 0.969).(b) Simulated trajectories of wage share (x 1 ) and employment rate (x 2 ) for the system in (a).
Subfigs.(a) and (c) in Fig 2 show a realization of the dynamics of the Keen system when the initial condition corresponds to a point in the region R d , and a finite equilibrium is reached.Subfigs.(b) and (d) show realizations when the initial condition belongs to R u .Under the latter setting, the system diverges as wages and employment go to zero while debt increases to infinity.Fig 3 shows regions R d (in gray) and R u (in white) for different levels of debt share (x 3 ).Figure (a) portrays the case for x 3 = 0.3, and figures (b), (c), and (d) present the cases for x 3 = 1, x 3 = 2.5 and x 3 = 5 respectively.

Fig 2 .Proposition 2 . 2 .
Fig 2. Keen model.Comparison of two Keen models under different initial conditions in the first and second rows of the panel.Subfigs.(a) and (b) portray the convergence to a stable finite equilibrium as the initial condition is given by x 0 = (0.75, 0.75, 0.1) > .Subfigs.(c) and (d) show the convergence to the undesired equilibrium, where the debt share grows infinitely, as the initial condition is given by x 0 = (0.8, 0.75, 0.1) > .https://doi.org/10.1371/journal.pone.0295859.g002

Fig 3 .
Fig 3. Finite and non-finite equilibria in the Keen model.Each sub-figure shows different combinations of initial conditions for x 1 and x 2 , when x 3 is fixed at some level.If the solution of the Keen model a particular combination of initial conditions reaches a desired finite equilibrium, then the cell is colored gray.In any other case, the cell remains white.The fixed levels of x 3 at each sub-figure are: (a) 0.03 (b) 0.1, (c) 2.5 and (d) 5. https://doi.org/10.1371/journal.pone.0295859.g003 Fig 4 illustrates the oracle reference path with its starting and ending points.

Fig 4 .
Fig 4. Gradient-based planner.Illustration of the Oracle reference path, generated using Algorithm 1.The path starts at x 0 = (0.5, 0.5, 0.8) > 2 R u , which is well outside the stable region of the Keen model, and moves towards x τ = (0.87, 0.87, 1.4) > 2 R d .The vector field of the Keen model is represented by the gray arrows in the figure to illustrate the dynamics of the system.https://doi.org/10.1371/journal.pone.0295859.g004 ffl ffl ffl ffl ffl ffl ffl {zffl ffl ffl ffl ffl ffl ffl ffl ffl ffl } ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl {zffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl } hðxÞ ; ð9bÞ

Fig 5 .
Fig 5. Block diagram representation of a control problem.The policy value u is the result of some designed policy generator, which takes as input a correction signal given by the difference between reference path ℓ and output y. https://doi.org/10.1371/journal.pone.0295859.g005 is the version of the Keen model that is utilized hereafter to apply standard linear control techniques through changes over input v. See Fig 6 for a detailed scheme of the control process.

Fig 6 .
Fig 6.Block diagram representation of the control problem for the Keen model.The diagram includes the specific components of the policy generator articulated via the designed control policy u: one-step-optimal-control and feedback linearization.The rest of the diagram is exactly as described in Fig 5.https://doi.org/10.1371/journal.pone.0295859.g006

Fig 7 ,
where the control starts at x 0 = (0.5, 0.5, 0.8) > and finishes in the vicinity of x τ = (0.87, 0.87, 1.4) > 2 R d .The same figure shows the vector field of the Keen model as the collection of gray arrows, portraying the dynamics of the system.As expected, when low latency is presumed, the tracking of the reference path is performed with barely any misalignments, as shown in panel (a) of Fig 7.However, if high latency is considered, as it is a more realistic scenario, tracking of the reference path is slow and achieved with a certain difficulty.In fact, Fig 7 differentiates the final point (x 0 t ) from the target point (x τ ) to

Fig 7 .
Fig 7. Control of Keen's model under two latency scenarios.Scenario (a) presumes low latency (τ 1 = τ 2 = 1000), while scenario (b) considers high latency (τ 1 = τ 2 = 1) in the control policies.As can be seen, in scenario (a), the target point (black) is almost on top of the final point (gray).In fact, the reference path (continuous black line) and the controlled path (gray dashed line) are very close to one another, indicating that the system can be controlled with ease.In scenario (b) we observe that the final (black) and target (gray) points are further from each other and that the controlled path does not follow the reference path with the same proximity as in the previous scenario, which resulted in a slower response to the applied control policies due to the added latency.Nonetheless, in both scenarios, we see that the Keen model is effectively controlled from an initial point (x 0 ) towards a small enough neighborhood of final destination (x τ 0 ) within the stable region of the Keen system.Once the control policy is released in such neighborhood the uncontrolled path follows a stable downward (continuous gray line).

Fig 8 .
Fig 8.Control of Keen's model under different latency scenarios.(a) shows the tracking of trajectories x 4 , x 5 and x 6 with low latency (τ 1 = τ 2 = 1000), and (b) shows the same for the case of high latency (τ 1 = τ 2 = 1).As can be seen, in the low latency case, the control policies are able to fluctuate more freely than in the high latency case, as the effect of the latter is the reduction in the velocity of the policies.In both cases, however, is noticeable the reduction in the rate of growth in labor productivity, as well as the increment in the rate of growth of the labor force, while maintaining positive values for the interest rates, as enforced by the conditions in equation(20).
appear explicitly.The functions hðxÞ; L f hðxÞ; . . .; L pÀ 1 f hðxÞ define a key local coordinate transformation for the system about x 0 and thus define a new set of coordinates z = T(x) where T : R n !R p , explicitly given by

Table 2 .
Initial conditions. https://doi.org/10.1371/journal.pone.0295859.t002 ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl fflffl {zffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl ffl } (35) is the case, one can defineu ¼ D À 1 ðv À L p f hðx 0 ÞÞð36Þand the linear system(35)can be re-written as: